Microwave Surface Impedance Measurements of the Electronic State and Dissipation 
of Magnetic Vortices in Superconducting Iron-Based LiFeAs Single Crystals 
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LiFeAs is one of the iron-based superconductors having multiple gaps with the possible sign 
reversal. To clarify how those novel natures affect the energy dissipation of magnetic vortices, we 
investigated the microwave surface impedance of LiFeAs single crystals under finite magnetic fields. 
The flux-flow resistivity enhanced rapidly at low magnetic fields, which is similar to the case of MgB2 . 
This is probably the consequence of the multiple-gap nature and the gap anisotoropy. This suggest 
that the sign-reversal is not important for the flux-flow even for multiple-gap superconductors. As 
for the electronic state, the vortex core of LiFeAs turned out to be "moderately clean" . Furthermore, 
the mean free path inside the vortex core was much shorter than that outside, and was close to the 
core radius. These results strongly suggest a process speciflc to the core boundary is important for 
a scattering mechanism inside the vortex core. 

PACS numbers: 74.25.Wx, 74.15.Lh, 74.25. nn 



Since the discovery of LaFcAs0i_2;Fa; with Tc = 26 K 
iron-based superconductors (SCs) have attracted lots 
of attention. Because multiple bands contribute to the 
Fermi surfaces and the magnetic phase exists in the vicin- 
ity of the superconductive phase in the phase diagram, 
it is expected that the mechanism of superconductivity 
of iron-based SCs is different from that of conventional 
SCs. New possibilities of superconducting gap structures 
based on the interband scattering, such as s^-wave @, 0| 
and s+^-wave [1, were suggested theoretically. Ex- 
perimentally, although this issue is under a debate Q, 
phase-sensitive experiments [tI-PI suggested that s^-state 
was realized in some materials of iron-based SCs. It is of 
great interest what the electronic structure and dynamic 
properties of vortices in such novel class of SCs are. 

As for conventional SCs, the quasiparticle (QP) exci- 
tation inside the vortex core has quantized energy levels 
with the spacing, AE ~ A^/E'p = Tiloq, where A and Ep 
are the size of the superconducting gap and the Fermi 
energy, respectively, and with those width, SE ^ h/rcoi-c, 
where Tcore is the relaxation time of QPs inside the vor- 
The ratio of these two energy scales, 
is a barometer of the quantum nature 
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of the electronic state inside the vortex core. Depending 
on this number, we have three regimes as i) the dirty 
regime (woTcoic ^ 1), ii) the moderately clean regime 
(luqTcoyc 1) and iii) the superclean regime (wq Tcore ^ !)■ 
It should be noted that Tcore is connected to the vis- 
cous drag coefHcient, ry, and the carrier density, n, as 
WoTcorc = ri/n-Kh [12}. 

According to Kopnin and Volovik (KV) [l^, the flux- 
fiow resistivity of a single-gap SC, p/, behaves in mag- 
netic fields, B, as 



Pi 

Pn 



A2 



B 



[B « B,2) 



(1) 



where pn, Bc2, Aq and (A^(0))fs are the resistivity in 
the normal state, the upper critical field, the maximum 
size of the superconducting gap and the angular average 
of the superconducting gap on the Fermi surface, respec- 
tively. This suggests that i) p/ in low B region increases 
linearly with B and ii) the gradient, a = Aq/(A^(0))fsj 
becomes larger than unity when A(9) depends on the 
angle 9. In fact, for an isotropic gap case, the Bardeen- 
Stephen (BS) theory [14] obviously obeys Eq.([T]). On the 
other hand, in nodal and modulated gap case, an en- 
hancement with a > 1 at low B region has been also 
observed experimentally This also suggests that 

the so-called "Volovik effect" (the effect of the Doppler 
shift on QPs disperse caused by the circulating supercur- 
rents) is not important for the flux-flow in low B region, 
although it succeeded to explain B dependences of the 
specific heat and the thermal conductivity. As for the 
2-band s++-wave SCs, such as MgB2 and Y2C3, a rapid 
enhancement of Pf{B) was observed This can 

be interpreted as the superposition of two linear B de- 
pendences corresponding to two bands 21|. Thus, pf{B) 
reflects the superconducting gap structure and its sym- 
metry. Therefore, it is very interesting how the flux-flow 
resistivity of the novel class of SCs behaves as a func- 
tion of B. However, the flux-flow of such novel SCs has 
not been investigated at all both theoretically and exper- 
imentally. Thus, it is a great challenge to investigate the 
flux-flow of iron-based SCs. 

We focus on a 111 material, LiFeAs. This material has 
the highest Tc of 18 K .22] among stoichiometric iron- 
based SCs, and single crystals with high quality (resid- 
ual resistivity ratio (i?i?i?)~50) can be obtained. The 
band calculation [i^l suggested that Fermi surfaces con- 
sist of two hole-like and two electron-like pockets around 
F-points and M-points, respectively. Nodeless multi- 
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pie superconducting gaps were observed by an angle- 
resolved photoemission spectrosc opy (ARPES) 
and a specific heat measurement [2^, superfluid-density 



data [27|, |28| showed that LiFeAs has nodeless multiple- 
gap structure. In addition to the phase sensitive exper- 
iment in Li-111 0, the electrical conductivity, ai [28t . 
estimated from the microwave surface impedance and 
the nuclear spin-lattice relaxation rate, 1/Ti do not 
show the so-called "coherence peak" below Tc. These 
strongly suggest that LiFeAs has the s^-wave gap struc- 
ture. Therefore, we can stand for the standpoint that 
Li-111 is an s±-SC. 

In this paper, we report the surface impedance of 
LiFeAs single crystals under finite magnetic fields, and 
discuss the electronic state inside the vortex core. It was 
clarified that the field dependence of the flux- flow of s^- 
state is similar to that of s+^-state, and that the vortex 
core of LiFeAs is moderately clean. The estimated mean 
free path of QPs inside the vortex core was found to be 
much shorter than that outside, and comparable to the 
core radius. This suggest that the mechanism character- 
istic to the core boundary plays an important role in the 
dissipative process inside the vortex core. 

LiFeAs single crystals were grown by a self-flux method 
[isl and were cleaved under Ar atmosphere in a glove box. 
Typical size of sample was 0.5x0.5x0.2 mm'^, and the de- 
magnetization coefficient estimated under ellipsoidal ap- 
proximation was about 0.58. These were of very high 
quality with RRR = /9dc(300 K)/pdc(Tc) ~ 45, and the 
dc resistivity behaved as pdc{T > Tc) ^ po + AT^ {po w 
30 ^ricm, A w 6.5 X 10~^ ^ficm/K^), which is typical 
of the Fermi liquid dominated by the electron-electron 
scattering. Since LiFeAs is moisture/atmosphere sensi- 
tive, samples were covered with Apiezon N grease during 
the measurement. We confirmed that Apiezon N grease 
does not affect results discussed below in a different com- 
parative experiment. 

The microwave surface impedance was measured by us- 
ing a cavity perturbation technique [s^ with a cylindrical 
oxygen- free Cu cavity resonator operated at u}/2tt ^ 19 
GHz in the TEqh mode. The Q-factor was Q > 6 x 10^, 
and the filling factor of samples was about 6 x 10^^. Both 
the external magnetic field up to 8 T and the microwave 
magnetic field were applied parallel to the c-axis. There- 
fore, we investigated the in-plane vortex motion. 

The surface impedance, Zg = Rs — iXg {Rs and Xg 
are the surface resistance and the surface reactance, re- 
spectively), is related to the resonant frequencies, ujg /2Tr 
and Wb/27r, and the Q- factors, Qs and Qb, as Rs = 
G(1/2Q, -l/2Qb), Xs = G{l-iJs/uJb) + C, where 
subscripts s and b represent the values measured with- 
and without the sample, respectively, and G, C are con- 
stants determined by the size and the shapes of the sam- 
ple and the resonator. The magnitudes of Rs and Xs 
are obtained by assuming the Hagen-Rubens relation, 
Rs = Xs = -y/ /iowpdc/2, in the normal state. 



Zs in the mixed state was calculated by Coffey and 
Clem (CC) [sij. Their calculation is based on the equa- 
tion of motion of the massless vortex, rjii + ku = ^qJ x 
z + f{t), where u is the displacement of a vortex, k is the 
pinning force constant, $o = h/2e — 2.07 x 10^^^ Wb is 
the flux quantum, J is the transport current density and 
z is the unit vector in the applied field direction. The 
effect of thermal fluctuations and the Hall effect are effec- 
tively included in random force, /(i), and rj, respectively 
for circulating microwave currents. At low temperature, 
the flux-creep contribution becomes negligibly small and 
the CC model leads to the relation 
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where Al is the London penetration depth, and Wcr/27r is 
the crossover frequency characterizing the crossover be- 
tween reactive- and resistive response, and s — p,oUjX^/pn 
which represents the normal-fluid contribution. One can 
assume that s to be negligible at low temperatures. Con- 
sequently, we obtain Wcr and pf from experimental data 
of Rs and Xg, by solving Eq.©. 

Figure [1] shows the magnetic field dependence of Zg at 
various temperatures. Good agreement between temper- 
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FIG. 1. (Color online) The magnetic field dependence of (a) 
tlie surface resistance, _Rs, and (b) the surface reactance, Xs, 
of a LiFeAs single crystal at 19 GHz up to 8 T at various tem- 
peratures. The curves and the open circles represent the data 
taken in the swept magnetic field (fixed temperatures) and in 
the swept temperature (fixed magnetic fields), respectively. 

ature swept data and magnetic-field swept data repre- 
sents that the magnetic field penetrates uniformly in the 
sample. With increasing magnetic field, both Rs and Xs 
increase monotonically. In particular, Rs shows a con- 
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vex upward behavior. We determine the zero-field super- 
conducting transition temperature, T°'^^°^ = 17 K, from 
the temperature dependence of Xg in zero magnetic field, 
which is in good agreement with the previously reported 
number in the same batch (ijj . 

The crossover frequency of Wci/Stt w 3 GHz obtained 
is larger than that of conventional SCs (w 100 MHz) [sl] 
but smaller than that of copper-oxide SCs by one order 
of magnitude 3^ 34 1. A similar value of Wcr has been 
reported in a 1111-type polycrystal (« 6 GHz) [sl]. The 
tendency that Wcr becomes small at high temperatures 
is consistent with a general description that the thermal 
fluctuation decreases the pinning force. 

Figure [2] shows the normalized flux-flow resistivity as 
a function of the normalized magnetic field. The flux- 
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FIG. 2. (Color online) The magnetic field dependence of the 
flux-flow resistivity pf{B) of the LiFeAs single crystal at sev- 
eral temperatures. The blue open circle is Pf{B) at T=1.8 
K obtained from temperature swept data. The gradient, a, 
expected in d-wave (with lines of node) SCs (a « 2) and in 
conventional s-wave SCs (a = 1) are also shown as dashed- 
and dotted lines, respectively. 

flow resistivity of LiFeAs single crystals increased lin- 
early with B, suggesting that the KV model is appro- 
priate even for this material. As for the gradient, a of 
LiFeAs is larger than that of the conventional s-wave case 
(a = 1) and smaller than that of the d-wave (with lines 
of node) case (a « 2). This enhancement of Pf{B) may 
be derived from one or both of two origins. First possible 
origin is based on the multiple-band nature. As for the 
2- bands SCs, such as MgB2 and Y2C3, the superposition 
of two linear dependences corresponding to two bands 
cause the flux-flow resistivity enhanced rapidly at low 
B [1^ [i^l- We can speculate that the 5-bands nature 
of LiFeAs probably induces the similar tendency. Sec- 
ond possibility is based on the gap anisotropy. Recent 
ARPES data suggests that some of the superconducting 
gaps have obvious 4- fold angle dependences [25|. Based 
on the KV model, this angle dependence of the super- 
conducting gap will make the gradient of Pf{B) larger 



than unity (a > 1). In any case, the magnetic field 
dependence of pf of LiFeAs are very similar to that of 
MgB2, implying that the s^-wave SC behaves essentially 
similarly to the s+^-wave SC so far as the flux-flow is 
concerned. The insensitivity of the fiux-fiow to the sign 
change for single-gap SCs has been already known for 
single-gap SCs; although the anisotropic s-wave SC and 
the d-wave SC differ from each other in the sign change 
of the order parameter, of both SCs show the B-linear 
dependence with a > 1 |16l4l8|. Our present result shows 



that the insensitivity shown up in the flux-flow is appli- 
cable also for multiple-gap SCs. 

Figure 13] shows the temperature dependence of the vis- 
cous drag coefficient, rj = ^qB/ pf. 77 is well fitted by 
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FIG. 3. (Color online) The temperature dependence of the 
viscous drag coeflicient, 77 = ^oB/pf. The green solid line is 
the expectation in the GL theory, ri{t) — ri{0){l — t'^)/{l + t^). 
The inset shows the temperature dependence of Bc2 obtained 
from from the temperature dependence of Ra (blue 

open circle), and that calculated from pn and rj (red open 
diamond). The solid line is eye guide. 

the expected temperature dependence in the Ginzburg- 
Landau (GL) theory, ri{T) = ri{0)[l - {T/Tc)^]/[l + 
{T/T^f]. From the fitting, we obtain 77(0) = (1.5±0.2) x 
10"^ Ns/m^. We can estimate the upper critical field 
as B,2{T) = p„(r)rKT)/$o, where p„(T) = po + AT^ 
is extrapolated pdc{T > Tc) to the temperature regions 
T < Tc- The result is shown in the inset. We obtain 
-Bc2(0) = 22 ± 4 T. Similar numbers were reported pre- 
viously in the same material 3B-41|. Considering the 
moisture/atmosphere sensitive nature of LiFeAs and the 
difference of RRR values among these crystals, we con- 
sider that it is within the range of individual differences. 

In figurcUl we discuss the relaxation time and the mean 
free path (mfp) of QPs inside the vortex core. Since 

'mo = 



Al (T) — pQe^ns{T)/m*, using the value m* 
5.2 — 6.3 (mo is the free-electron mass) [i^, |43| and 
Al(0) = Xs{0)/fJ,oLd ^ 390 nm, we estimate the carrier 
density n « ns{0) = (9.6 — 11.7) x lO^'^ cm~^, which gives 
i^o''coiG = 0.4 — 0.5. This shows that the vortex core of 
LiFeAs is in the moderately clean regime. Furthermore, 
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by using the number — ?iwo/2 = —0.9 meV observed 
in a recent scanning tunneling microscopy/spectroscopy 
(STM/STS) study "44], we obtain the relaxation time of 
QPs inside the vortex core, rcorc(l-8 K) « 0.15 ps. This 
value is quite different from that outside (w 10 ps) [ist . 
and is even smaller than that in the normal state (~ 0.6 
ps). These are shown in Figure |4lJ a). From the relaxation 
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FIG. 4. (Color online) (a) Temperature dependences of re- 
laxation times and mean free paths of QPs in several states. 
Symbols are those outside the vortex core (blue open circle), 
in the normal state (green open triangle) and inside the vortex 
core (red solid circle), respectively, (b) Temperature depen- 
dences of the mfp inside the vortex core (red circle) and the 
coherence length calculated from Bc2{T) (black triangle). 

time, we found that the mfp of QPs inside the vortex core 
to be /core — VFTcore ~ 40 A, where vp ~ 2.6 x 10^ m/s 
is the Fermi velocity, which is estimated from STM/STS 
13] and ARPES data [13, Hi] . Again, this value is much 
shorter than that outside the core, ^Meissner- In particu- 
lar, as shown in Figure I^Jb), ^corc is comparable to the 
coherence length, ^, estimated from Bc2- We checked 
the repeatability in another single crystal of LiFeAs, and 
the results were consistent with those described above. 
In addition, we performed the same measurements in 
LiFe(As,P) single crystals, which was at most 3 % P- 
substituted, and we obtained the similar results. 

The short mfp of QPs inside the vortex core 
was also observed in many copper-oxide SCs, such 
as YBa2Cu307_a;, Bi2Sr2CaCu20y and La2-a;Sra;Cu04 
IB" 47 1- In these cuprate, the mfp inside the vortex core 
is also much shorter than that outside and rather close 
to the core radius, Im ^ ^corc ~ Similarly, in Y2C3 
[20| . which is one of the 2-gap SCs with isotropic s-wave, 
the mfp inside the vortex core is limited to the coherence 
length, Zcorc ^ It is surprising that similar tendency 
was observed among many different SCs with different 
gap structures, pairing mechanisms and electronic struc- 
tures. Since the relation, Zcorc ^ was obtained, one 
can consider that a scattering process which is specific 
to the core boundary contributes to the additional dissi- 
pation in the vortex core as was originally considered by 



core boundary is crucial even in the flux-flow of mod- 
erately clean SCs, and theoretically showed that there 
is extra energy dissipation at low frequencies because of 
the presence of a collective mode. However, it is not yet 
clear whether this mechanism can explain the large dis- 
sipation observed in our experiments quantitatively at 
present. Systematic study of the frequency dependence 
of the in-core dissipation will clarify the validity of Es- 
chrig's model. On the other hand, according to Tinkham 
[Hot Nozieres-Vinen- Warren [3 lEH , the relaxation 
time Tgap = h/ Aq which is characteristic of the moving 
vortex, has been considered. For LiFeAs, Tgap = 0.2 ps is 
comparable to obtained Tcorc- In order to clarify the va- 
lidity of these models, studies of the gap-size dependence 
of Tcorc is needed. 

In conclusion, we investigated the microwave surface 
impedance of LiFeAs single crystals under finite magnetic 
fields. The magnetic field dependence of the flux-flow re- 
sistivity of new class of superconductors having multiple 
gaps with the possible sign reversal became clear. The 
flux-flow resistivity increased linearly with the magnetic 
field, as was suggested by Kopnin-Volovik. Particularly, 
the gradient at low fields was larger (smaller) than that 
of conventional s-wave superconductors (d-wave super- 
conductors with lines of node). This is probably the 
consequence of the multiple-gap nature and/or the gap 
anisotropy. This also suggests that the flux-flow resis- 
tivity is insensitive to the sign reversal of the order pa- 
rameter on different Fermi surfaces. As for the electronic 
state, the vortex core of LiFeAs was estimated to be the 
moderately clean. The mean free path of quasiparticles 
inside the vortex core was much shorter than that out- 
side, and comparable to the core radius, suggesting the 
importance of the Andreev reflection at the core bound- 
ary. Such a tendency was observed also in many other 
superconductors, and systematic studies will clarify the 
dissipative mechanism inside the vortex core. 

We thank Tetsuo Hanaguri for showing us many un- 
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thank Masashi Takigawa for providing us LiFeAs single 
crystals, and Yusuke Kato for valuable comments. 



Nozieres and Vinen for clean SCs [48[. Indeed, Eschrig 
et al. |4£] discussed that the Andreev reflection at the 
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